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We present the complete zero temperature phase diagram of a model for ultrathin films with perpendicular 
anisotropy. The whole parameter space of relevant coupling constants is studied in first order anisotropy approx- 
imation. Because the ground state is known to be formed by perpendicular stripes separated by Bloch walls, a 
standard variational approach is used, complemented with specially designed Monte Carlo simulations. We can 
distinguish four regimes according to the different nature of striped domains: a high anisotropy Ising regime 
with sharp domain walls, a saturated stripe regime with thicker walls inside which an in-plane component of 
the magnetization develops, a narrow canted-like regime, characterized by a sinusoidal variation of both the 
in-plane and the out of plane magnetization components, which upon further decrease of the anisotropy leads to 
an in-plane ferromagnetic state via a spin reorientation transition (SRT). The nature of domains and walls are 
described in some detail together with the variation of domain width with anisotropy, for any value of exchange 
and dipolar interactions. Our results, although strictly valid at T = 0, can be valuable for interpreting data on 
the evolution of domain width at finite temperature, a still largely open problem. 
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I. INTRODUCTION 

The magnetic phases of ferromagnetic thin films with per- 
pendicular anisotropy have been the subject of intense exper- 
imental 1-7 , theoretical 8-14 and numerical 15-20 work in the last 
20 years. Magnetic order in ultrathin ferromagnetic films is 
very complex due to the competition between several different 
energy contributions, the most prominent being exchange and 



dipolar interactions, together with a strong influence of shape 
and magnetocrystaUine anisotropies of the sample. These 
in turn are very susceptible to the growth conditions of the 
films^i. 

A widely used model that contains the main ingredients of 
ultrathin film magnetism is the 2D dimensionless Heisenberg 
Hamiltonian: 
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where Si are classical unit vectors, the exchange and 
anisotropy constants are normalized relative to the dipolar 
coupling constant (S = J/ft./rj = K/Vt), < i,j > stands for 
a sum over nearest neighbors pairs of sites in a square lattice, 
(i, j) stands for a sum over all distinct pairs and r-jj = \ fi — fj \ 
is the distance between spins i and j. 

At low temperatures and strong enough perpendicular 
anisotropy, the presence of a striped phase (i.e., a modulated 
pattern of local perpendicular magnetization with a well de- 
fined stripe width h) is well established and is the ground state 
of the syste m 9 ' 10,22 . In the limit of strong uniaxial anisotropy 
domain walls are sharp and the energy cost for deforming 
or moving a domain wall is large. Nevertheless, even when 
the mechanism by which the width of domains adjusts is not 
well understood, the stripe width varies with the effective 



anisotropy. When the thickness of the films (or the tempera- 
ture) grows, the effective perpendicular anisotropy is reduced 
in films of a few monolayers, and magnetostatic energy be- 
comes important, inducing the magnetization to develop an in- 
plane component. Domains become narrower, walls become 
wider and are of Bloch type at low temperatures 4 ^, until the 
system goes through a Spin Reorientation Transition (SRT) 
when anisotropy and dipolar energies cancel Around 
the SRT line a canted state may develop, where the magneti- 
zation presents a finite in-plane component together with the 
perpendicular modulation. The extension of the canted state in 
parameter space strongly depends on the nature of the relevant 
anisotropies. For some systems, like Co/Au(lll) , it seems 
necessary to go beyond the first order anisotropy approxima- 
tion of the model (Q~|). A second order anisotropy energy is 
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responsible for a canted state in a large portion of the phase 
diagra m 14 ' 16 ' 25 . For other systems, like Fe/Cu(001), the first 
order anisotropy seems to be enough to describe the relevant 
physic s 27 ' 28 . In this case, as described by the model (Q~|i, the 
canted state is restricted to a narrow region of parameter space 
around the SRT, as expected from general considerations 16 , 
and reported in simulations at finite temperature 29 . Finally, 
when the dipolar anisotropy exceeds the magnetocrystalline 
one, the system enters an in-plane ferromagnetic state. 

In this work we extend previous calculations 9 ' 13 ' 30 and com- 
pute the complete phase diagram in the (5, 77) space of Hamil- 
tonian ((TJ, at T = 0. We also improve upon previous re- 
sults by considering different kinds of domain walls (sinu- 
soidal, hyperbolic tangent, sharp walls), as appropriate for 
each regime in parameter space. We consider only straight 
domains, (domains in which the spin orientation can be modu- 
lated along the x direction but is constant in the perpendicular 
direction y) separated by Bloch walls, i.e., walls in which the 
magnetization stays inside the yz plane. The local magnetiza- 
tion vector inside the domains may be tilted at an angle 9 with 
respect to the plane normal (z axis). Within these assumptions, 
we obtain the complete phase diagram, the variation of the an- 
gle 9 and the behavior of the width of domains and walls in 
the whole parameter space (6, 77). This allows, e.g., to obtain 



the crossover between Heisenberg (extended walls) and Ising 
(sharp walls) regimes. We verify that domain width adjust- 
ment with varying anisotropy is only possible in the Heisen- 
berg regime, domain width being fixed in the Ising regime 
for any value of the parameters (<5, 77). Analytical calculations 
are complemented with Monte Carlo simulations specially de- 
signed for the present purposes, as explained in the appendix. 



H. ZERO TEMPERATURE PHASE DIAGRAM 

We consider a square lattice with N = L x L sites, charac- 
terized by the integer indexes (x,y), where —L/2 < x < L/2 
and —L/2 < y < L/2, in the limit L — ► 00. Hence, the in- 
dex i in Eq.(Q~|i denotes a pair of coordinates (x, y). We con- 
sider only uniformly magnetized solutions along every verti- 
cal line of sites, i.e. or x ,y\ = M(x), Vy and allow only Bloch 
walls between domains of perpendicular magnetization, i.e. 
M x (x) = Vx. Yafet and Gyrogy (YG) showed that for these 
types of spin configurations the energy per spin can be mapped 
onto the energy of a one dimensional XY model 9 . The energy 
difference between an arbitrary magnetization profile M(x) 
and a uniformly in-plane magnetized state is then given by: 
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where 5' = S - 2 a, k' = 77 - 3g, c x = 0.01243 . . ., c 2 = 
0.07276 . . ., g = 1.202057 ... and 

C = C [M y (x)} = 2(c 2 - c x )- M y {x) M y (x + 1) (3) 

X 

Although small, this correction term makes a non negligi- 
ble contribution when the domain walls are of the same or- 
der of the lattice constant. This happens for small values 
of 5 (S < 5), where both the stripe and wall widths are of 
the order of a few lattice spacings. For larger values of 5 
it is reasonable to assume a smooth magnetization profile 9 



M v (x + 1) w M v (x), so that the correction © can be ab- 
sorbed into the anisotropy term in Eq.(f2]i, replacing n' — > n = 
r]-3g + 2(c 2 - ci). 

Now consider a stripe-like periodic structure of the magne- 
tization profile with period 2h, M z (x + h) = —M z (x). Using 
a Fourier expansion: 



M z (x)=M J2 &™cos(^p), (4) 

m=l,3,... 

the energy (O can be written as 9 
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where 4>(x) is the angle between M(x) and the z axis and 
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Now we look for the minimum of Eq.© for different val- 
ues of 5, r\. We propose different striped magnetization pro- 
files M z (x) and compare the energies obtained by minimizing 
Eq.© for each profile with respect to variational parameters. 

We first consider a profile as proposed by YG, that is con- 
stant | M z (x) = Mq inside the stripes with a sinusoidal varia- 
tion inside the walls between stripes (see Fig.l in Ref. 9 ). This 
will be called "sinusoidal wall profile approximation" (SWP). 
In order to allow for canted profiles, we take Mo = cos 9, 
where 9 is the canting angle, i.e. we define it as the mini- 
mum angle of the local magnetization with respect to the z 
axis. In Ref. 9 this variational problem was solved for Mq = 1 
in the continuum limit, i.e. when h ^> 1 and the wall width 
w > l,so that the profile can be considered a smooth function 
of x. While this approximation is expected to work well for 
large enough values of 8, it breaks down for relatively small 
values of it, where the discrete character of the lattice has to 
be taken into account. However, the variational problem for 
that range of values of 6 can be solved exactly (although nu- 
merically) by minimizing Eq.© with respect to the integer 
variational parameters h and w and continuous parameter 9. 
In other words, for every pair of values (8, rf) we evaluate the 
energy Eq.© for the sinusoidal profile with different combi- 
nations of h = 1, 2, . . . and w — 1,2, . . . within a limited set. 
For every pair of values h, w,we look for the value of 9 that 
minimizes the energy with a resolution A9 = 0.01 and com- 
pare all those energies. This calculation is feasible for values 
up to 5 = 10, for which the maximum value of h (bounded by 
the stripe width in the 77 — > 00 limit) remains relatively small 
(smaller than h = 140). Some results for <5 = 15 close to the 
SRT were also obtained. All the results of this calculation are 
compared against Monte Carlo (MC) simulations. Details of 
the MC method used are given in AppendixlAl Through these 
calculations we obtain a zero temperature phase diagram for 
low values of 8. 

Before presenting the results, it is important to introduce 
some notations and definitions of different types of solutions. 
We distinguish between four types of solutions. If the mini- 
mum energy solution corresponds to w = 1 and 9 = (within 
the resolution A9), we call this a Striped Ising Profile (SIP), 
i.e. a square wave like profile. If 9 = but w > 1, we call 
this a Saturated State (SS). These states only show a finite 
in-plane component of the magnetization inside the walls. If 
< 9 < 7r/2 the solution is a canted-like state. Finally, if 
9 = ir/2 (Mo = 0) we have a Planar Ferromagnet (PF). 

The zero temperature phase diagram for small values of 5 
(8 < 5) is shown in Fig[T] For relatively large values of 77 the 
minimum energy configuration is always the Ising one (SIP), 
with a stripe width independent of 77. For small values of r\ 
the minimum energy configuration is the PF, with a spin re- 
orientation transition line (SRT), either to the Ising state for 
h < 3 (8 ~ 2)) or to a canted-like one for h > 3 (8 > 2). No 
Saturated State configurations are observed for 8 < 6. 

Inside the canted region, a strong stripe width variation with 
the anisotropy is observed at constant 8 . Note that the ver- 
tical lines that separate Ising striped states with consecutive 
values of h bend inside the canted region and become almost 
horizontal as 8 increases. Hence, the exponential increase 




FIG. 1: (Color online) Zero temperature phase diagram for small 
values of 8. Black filled symbols and black solid lines: MC simula- 
tions. Open red symbols: SWP approximation. Squares and contin- 
uous black lines correspond to transition lines between striped states 
of different width. The shaded region corresponds to a canted-like 
state (0 < 6 < n/2). Triangles are transition lines between Pla- 
nar Ferromagnet (6 = n/2) and SMCP States (Spin Reorientation 
Transition line). Circles mark transitions between the canted-like 
and the Striped Ising state (6 = and w = 1). Notice the excellent 
agreement between the MC and SWP calculations close to the SRT, 
while the SWP approach underestimates the transition line between 
the canted-like and Ising Stripes states. The dotted line corresponds 
to the continuum approximation of YG for the SRT (Eq.Q). 



of h with S in the Ising region (vertical lines) changes to an 
exponential increase with r\ inside the canted region (curved 
lines on the right of Fig JTJ- It is important to note here that 
the canted region in this system corresponds almost every- 
where (except close to the crossover to an Ising striped state) 
to a regime in which the stripe width and walls are of the al- 
most equal, which means a pure sinusoidal magnetization pro- 
file. In this sense it has a different character than the "true" 
canted phases obtained in systems with non-zero higher or- 
der anisotropics 14,16 , where well defined domains show a fi- 
nite in-plane magnetization component. In the present case 
the canted like states are characterized by a sinusoidal vari- 
ation (with wave length 2h) in both the in-plane and the out 
of plane magnetization components, without well defined do- 
mains (see an example in Fig©. Hence, there are not truly 
"stripes", but a sinusoidal modulated state or Single Mode 
Canted Profile (SMCP). 

We also find an excellent agreement between the sinusoidal 
wall approximation (or SMCP) and the MC results, except 
close to the transition between the Ising and the canted-like 
states. Such disagreement is due to the fact that the actual 
wall is not well described by a sinusoidal profile far away of 
the SRT line, as will be shown later. 

For large enough values of 8 the variational problem for the 
SWP can be solved in a continuum approximation introduced 
by YG 9 , giving a set of coupled non-linear equations for the 
stripe width h, the ratio between the stripe and the wall widths 
A = w/h and the canted angle 9. In the limit A — > 1 those 
equations can be solved analytically predicting a SRT at the 
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FIG. 2: (Color online) Canted angle (circles) and in-plane magneti- 
zation (squares) as a function of r\ for S — 7.5. Filled black symbols 
correspond to MC calculation. Open red symbols corresponds to the 
discrete SWP approximation, while the red doted and dashed lines 
correspond to the continuum (YG) approximation of the SWP. Con- 
tinuous black and red lines are only a guide to the eye. 
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Indeed, from our MC simulations, we observe that the 
range of values for which canted-like states have the mini- 
mum energy gradually shrinks as 6 increases, being replaced 
by a saturated state for values of 77 above certain threshold. 
This can be observed in Figj2] where we show the behavior of 
the canted angle and the in-plane magnetization component 
M|| = (1/L) J2 X My{x) as a function of 77 for S = 7.5. The 
Monte Carlo data shows the existence of a wide range of val- 
ues of 77 for which the canted angle is zero while M | 7^ 0, 
meaning that the non null in-plane components are concen- 
trated inside the walls. In other words, in that region we have 
a saturated state with thick walls w > 1. Notice also that 
the SWP approach completely fails to describe those states. 
Moreover, we observe from our MC simulations that the SWP 
cease to be the minimum energy solution for values of 77 rela- 
tively close to the SRT, before the saturated state sets up (see 
Fig|2]i. This effect becomes more marked as 5 increases. 

The departure of the magnetization profile from the SWP 
for large values of 77 and 5 is expected from micromagnetic 
theory, which in that limit predicts that the wall structure will 
be dominated by the interplay between anisotropy and ex- 
change, leading to an hyperbolic tangent shape of the wall 23 . 
Hence we considered a periodic magnetization profile with 
hyperbolic tangent walls (HWP) defined, for a wall centered 
at x = 0, by 



Vsrt(S) = a - — 



(7) 



with a = 7r 2 /3 + 3g — 2(c2 — ci). The line Eq.© is also 
depicted in FigQ] Notice the disagreement between the con- 
tinuum approximation and the exact one for 6 < 5. This dis- 
crepancy becomes smaller than 1% only for 5 > 7. 

For arbitrary values of 77 and 6 the equations for h, A and 8 
can be solved numerically. From the numerical solutions we 
found that the range of values of the anisotropy 77 for which 
the canted angle is appreciable different from zero within this 
approximation is strongly depressed as S increases. For val- 
ues S ~ 100 the canted-like configurations almost disappear, 
except very close to the reorientation line, as already reported 
by Politic. 



M z {x) = M tanh 



for -h/2<x< h/2, (8) 



where Mo = cos 8 as before. In the large S limit, assuming a 
smooth profile h ^> 1 and l w 1, the anisotropy energy can 
be expressed as: 
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The exchange energy can be obtained in a similar way: 
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The dipolar energy can be calculated using Eq.©. The 
Fourier coefficients for the profile (0 can be computed using 
the approximation 



tanh(x) 



(1- 
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(11) 

This leads to an expression for the total energy as a function 
of the variational parameters h, 8 and l w that can be mini- 



mized numerically. Comparing the minimum energy solution 
for the SWP and the HWP we obtain the crossover line be- 
tween sinusoidal and hyperbolic wall structure shown in Fig[3] 
(dashed line). Above that line the HWP has always less energy 
than the SWP. We also calculated the transition line between 
the canted-like and the saturated states by setting the condi- 
tion 8 = 0.01, to be consistent with the criterium used in 
the MC calculations. The results are shown in Fig [3] together 
with the SRT line Eq.©, and compared with MC calculations 
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FIG. 3: (Color online) Zero temperature phase diagram for large val- 
ues of 8. The shaded region corresponds to the canted-like states. 
Symbols correspond to MC simulations and lines to theoretical re- 
sults. The dashed line correspond to the crossover between sinu- 
soidal and hyperbolic wall structure. The lower line (blue) corre- 
sponds to Eq.Q. The middle line (red) is obtained from the HWP 
minimum energy solution with 8 = 0.01. The upper line (black) 
corresponds to Eq.l|19t. Typical magnetization profiles obtained by 
MC are shown for every region of the phase diagram. 



with 7 = 7r 2 /3 — k. Minimizing Eq.dTZb with respect to the 
variational parameters h and l w leads to: 
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3^ 
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with 
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(14) 



in agreement with a derivation made by Politic. 

With the previous calculation we can also estimate the tran- 
sition line between the saturated and the Ising Striped state. In 
the large h limit the energy for a SIP, i.e. for 



, , , / if 0<x<h/2 

= {ir if h/2<x<h (15) 

the energy can be easily calculated from Eq.©. The Fourier 
coefficients are: 



up to 5 = 15. The excellent agreement with the MC results 
gives support to the analytic approximations. Only between 
the dashed and the continuous lines we found truly canted 
states (i.e., states with h <C w), although they parallel com- 
ponent is rather small (9 < 0.1). In this sense, the region 
enclosed by both lines marks a crossover between the SMCP 
and the Saturated states: as r\ increases domains emerge grad- 
ually, the walls change from sinusoidal to hyperbolic shape 
and the canted angle goes to zero. 

For large values of r\ the exponential increase of h makes it 
cumbersome to apply the previous approximation for the cal- 
culation of the dipolar energy. Instead of that, we can use the 
following heuristic argument to obtain a reasonable approxi- 
mation. The main error introduced by the SWP approach is 
in the exchange and anisotropy contributions to the energy. 
Since the main contribution to the dipolar energy is given 
by the interaction between domains, we can assume that the 
dipolar contribution of the wall is relatively independent of 
its shape. Hence, we approximate the dipolar contribution by 
the SWP expression obtained by YG 9 taking w = f l w (/ is 
a fitting parameter of order one to be fixed later) in the limit 
A <C 1 (l w /h -C 1). We compare the energy obtained within 
this approximation with that obtained using the Eq.(fTTT> for 
different values of the system parameters. We verified that the 
error made by taking / = 4 is always smaller than 1% for 
h/l w > 20. We also observe that the best agreement with the 
MC results is obtained for / = 4. Assuming then Mq = 1, 
the total energy per spin (relative to the parallel magnetized 
state) for the HWP can then be approached by 
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Using Eq.© the dipolar energy is then given by 
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where /3 = j e + In 4 — 1, j e f=s 0.577216 is the Euler gamma 
constant and ip(x) is the digamma function 31 . The energy per 
spin respect to the in-plane magnetized state is then given by 
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Minimizing Eq.([T8l with respect to h leads to the equation 
S'/2 - (3 = F(h), where F(h) = ij)(h) - h ip'(h) ~ Inh - 1, 
thus recovering the known result h ~ e 5 / 2 . Comparing the 
energies, we find that the HWP has less energy than the Ising 
state for any value of rj. Eq.(IT3l shows that the stripe width 
variation in the Saturated state is determined by the change 
in the wall width as the anisotropy increases. Hence, h will 
change until the wall width reaches the atomic limit, i.e. for 



In 



1, where Eq.([T3l recovers the Ising behavior h 



J/2 



Imposing the condition l w = 1 to Eq.(fT4l we obtain the tran- 
sition line between the Saturated and the Ising Stripes states: 



1 7T 2 

V= ^S-2+ — + 3g-2(c 2 - Cl ) 



(19) 



which is also shown in Fig|3] in complete agreement with the 
MC results. 
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FIG. 4: (Color online) Comparison of the T — stripe width h vs. 
r\ obtained within the different methods for 5 = 10. Symbols corre- 
spond to MC simulations. Full black line corresponds to the asymp- 
totic approximation for the HWP given by Eqs.l|13t and The 
red dashed line corresponds to the variational solution of Eqs.l(9j and 
UOI l using the approximation jilt for the Fourier coefficients in the 
dipolar energy. The blue dotted line corresponds to the continuous 
SWP. The vertical dashed line corresponds to the transition between 
Saturated and Ising Stripes states given by Eq.l|19t. 



In Fig|4]we compare the equilibrium stripe width h as a 
function of 77 obtained within the different approximations 
used in this work for 5 = 10 and with the MC simulations. 
Notice that the asymptotic approximation for the HWP given 
by Eqs.dOll and ( fT4l i shows a better agreement with the MC 
results than using the approximation (fTTT i for the Fourier coef- 
ficients in the dipolar energy. This is because we adjusted the 
fitting parameter / to optimize the agreement with the MC re- 
sults at low values of i5. The discrepancy between both (hyper- 
bolic) approximations becomes negligible in the large 5 limit. 



III. DISCUSSION AND CONCLUSIONS 

The main results of this work are summarized in FigsQ]and 
|3] which display the complete zero temperature phase dia- 
gram of the model defined by the Hamiltonian ([TJ. Work- 
ing upon reasonable assumptions for the ground states, like 
perfectly straight modulations in one dimension and Bloch 
domain walls, we analyzed minimum energy configurations 
combining a variational analysis with Monte Carlo results. 
We found four qualitatively different kinds of solutions: a pla- 
nar ferromagnet for small anisotropics, a Single Mode Canted 
Profile (characterized by a sinusoidal variation of both com- 
ponent of the magnetization and varying wave length) close 
to the SRT and two types of perpendicular striped states for 
large values of the anisotropy: a saturated state in which the 
in-plane component is restricted to the domain walls, and an 
Ising stripe state with sharp walls for large anisotropics. 

The SMCP and saturated states give valuable information 
on the behavior of the stripe width (or the wave length in 
the SMCP case) as the anisotropy and exchange parameters 



change. We find that stripe width variation is directly associ- 
ated to the presence of finite width domain walls. For large 
enough values of the anisotropy 77 the ground state of the sys- 
tem is always an Ising Striped state, no matter the value of the 
exchange coupling 5. In those states domain walls are sharp, 
the stripe width is completely independent of 77 and grows ex- 
ponentially with the exchange coupling. 

At the SRT the system passes through canted-like states 
(mostly SMCP)as the anisotropy increases, although the range 
of values of 77 where the canted angle is different from zero 
narrows as S increases. For instance, the exchange to dipolar 
coupling ratio in fee Fe based ultrathin films can be roughly 
estimated to be S ~ 100 (considering a cubic bilayer of 
Fe/Cu(100), where 4 the exchange coupling Jp e ~ 30meV, 
the lattice constant dFe ~ 2ML and 32 /jp e ~ 3/is). For 
5 ~ 100 the anisotropy interval for the canted-like states is 
approximately A77 = 77 — t]srt ~ 0.2. 

For 6 < 6 the SMCP has the minimum energy in a rather 
extended region of the phase parameters space, close to the 
SRT. The wave length (or "stripe width") of those states 
presents a strong variation with the anisotropy, directly cor- 
related with an increasing canted angle. According to YG ap- 
proximation, SMCP solutions are expected only close to the 
SRT. We found that the magnetization profile maintains the 
sinusoidal shape as the anisotropy increases. Above certain 
value of 77 the wall profile changes to a hyperbolic tangent 
shape, while the magnetization inside the domains becomes 
fully saturated. 

For 5 > 6 the ground state is given by the Saturated State, 
except very close to the SRT. A similar effect (i.e. a crossover 
between a sinusoidal and a saturated magnetization profile) 
has been observed in room temperature grown fee Fe/Cu(100) 
ultrathin films, as the temperature decreases from T c , even 
though those systems do not present SRT—. 

In the Saturated state, the stripe width increase with 77 is 
directly related to the wall width decrease through the rela- 
tion h ~ e 5 / 2i ™. The wall width in turn is determined by 
the competition between exchange and anisotropy. Once the 
anisotropy is large enough that the wall width reaches the 
atomic limit, h growth stops. One may wonder whether a sim- 
ilar mechanism could be behind the stripe width variation with 
temperature, where a saturation is observed at low tempera- 
tures. Nevertheless, in this case other effects, like extremely 
slow relaxation can be responsible for the observed saturation. 

Besides its direct application to real systems, knowing the 
ground state of this system for arbitrary values of the exchange 
coupling is of fundamental importance to have a correct inter- 
pretation of Monte Carlo simulation results. Being one of the 
most powerful tools to analyze these kind of systems at the 
present (specially at finite temperatures), it is basically lim- 
ited by finite size restrictions, which implies relatively small 
values of S (the characteristic length h of the problem grows 
exponentially with 8 at low temperatures). 
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Appendix A: Zero Temperature Monte Carlo Technique for 
striped domain patterns 

In order to have an independent computation of the striped 
profiles which minimize the energy, we implemented Monte 
Carlo simulation with a simulated annealing protocol and 
Metropolis algorithm. To compare against the analytical so- 
lutions we looked for minimum energy magnetization profiles 
among those characterized by periodic straight domains with 
Bloch walls. Hence, the problem is basically one dimensional 
and we could restrict the search to one dimensional patterns 
over the x direction fixing Sf = Vi and imposing periodic 



boundary conditions (PBC) in both the x and y directions. In 
other words, we simulated a lattice with L x x L y with L y = 1 
and PBC, which are implemented by means of the Ewald sums 
technique. 

The temperature was then decreased down to very low tem- 
peratures at a constant rate T(i) = To — rt, where time is 
measured in Monte Carlo Steps, To is the initial temperature 
and r is the cooling rate. For all the range of parameters of 
this work, we choose Tq = 1 and r = 1CP 4 and the simulation 
protocol was repeated 100 times using different sequences of 
random numbers in order to minimize the possibility of trap- 
ping in local minima. The results were independent of the ini- 
tial spin configuration we choose at To. For every set of values 
of (6, if) we checked the results for different values of L x in 
order to avoid artificial frustration. We also performed some 
comparisons with MC results in a square L x = L y lattice with 
PBC using the same annealing protocol and the results were 
indistinguishable. This ansatz allowed us to obtain MC re- 
sults for values of 6 up to S = 10 (for which the maximum 
equilibrium value is h = 140). 
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